Quantum Phases of the Extended Bose-Hubbard Hamiltonian: Possibility of a 
Supersolid State of Cold Atoms in Optical Lattices 
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Cold atom optical lattices typically simulate zero-range Hubbard models. We discuss the theoret- 
ical possibility of using excited states of optical lattices to generate extended range Hubbard models. 
We find that bosons confined to higher bands of optical lattices allow for a rich phase diagram, 
including the supersolid phase. Using Gutzwiller, mean field theory we establish the parameter 
regime necessary to maintain metastable states generated by an extended Bose-Hubbard model. 
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Bose condensed cold atom systems in optical lattices 
are increasingly serving as beautiful (and practical) 
laboratories for studying quantum phases and quantum 
phase transitions in strongly correlated model Hamilto- 
nians of great intrinsic interest. In particular, the very 
low temperature; the absence of disorder, dirt and de- 
fects; and (essentially) complete control over the system 
parameters (and therefore the applicable Hamiltonian) 
combine to make cold atoms in optical lattices an ideal 
system to experimentally test the predictions of various 
interacting quantum Hamiltonians which originated as 
model (and often rather unrealistic) descriptions for con- 
densed matter physics problems. For example, the Mott 
insulating phase and the superfluid phase of a Bose Hub- 
bard model have been demonstrated in the cold atom 
optical lattices 0, Q ■ 

In spite of the impressive success of the cold atom sys- 
tems in studying the quantum phases of strongly corre- 
lated Hamiltonians, there has been one important limi- 
tation. Cold atoms in optical lattices usually represent 
essentially zero-range systems where the correlation (or, 
equivalently, the inter-particle interaction) is effectively 
on-site only, being parameterized by a single interaction 
energy U (the so-called Hubbard U), so that the sys- 
tem Hamiltonian is the Hubbard model characterized by 
a single dimensionless coupling parameter t/U, where t 
is the quantum tunneling or hopping strength. In this 
letter we propose a specific and practical scheme to gen- 
eralize the cold atom Hamiltonian to an extended Hub- 
bard model, where both on-site (U) and longer-range (V") 
inter-particle interactions compete with the kinetic en- 
ergy giving rise to a rich quantum phase diagram which 
should be experimentally accessible. We focus, in par- 
ticular, on bosonic systems though the method described 
here is applicable to spinful, fermionic systems as well. In 
the bosonic case, our proposed scheme may lead to den- 
sity wave and supersolid quantum phases in addition to 
the "usual" Mott insulating and superfluid phases. The 
key idea in our work, enabling the realization of an ex- 
tended Hubbard atomic system, is that one could, by 
utilizing proper laser excitations of individual cold atom 



states in the confining potential 0.0. use, in principle, 
the excited confined states (rather than the lowest level 
in each individual optical lattice potential minimum) to 
form the interacting system. Such a system would have 
a natural extended Hubbard description rather than an 
on-site description. We theoretically obtain the quantum 
phase diagram of such a system and predict the exciting 
possibility of coexisting density wave and superfluid or- 
der, i.e. the supersolid quantum phase, in the U — V 
parameter space of a realistic, extended Bose-Hubbard 
model. Our proposed system should also have important 
relevance to topological quantum computation |fj in op- 
tical lattices which has been shown to be feasible with 
cold atom lattices provided an extended range interac- 
tion (V 7^ Q) applies. 

Cold atomic gases confined to optical lattices offer the 
unique opportunity to directly probe novel states of mat- 
ter, including the supersolid. In comparison, experi- 
mental evidence for supersolid order in He now exists 
7], though conclusive identification using current exper- 
imental techniques remains elusive In bosonic cold 
atom optical lattices, coherencepeaks in multiple mat- 
ter wave interference patterns ||j at half the reciprocal 
lattice vector would provide strong evidence for super- 
solid order. Recent proposals suggest that nearby Mott 
(and therefore density wave) order may also be dire ctly 
observed, but through structure in noise correlations [lflj 
or through Bragg spectroscopy [Tl) . 

We begin with the following second quantized Hamil- 
tonian describing bosons in an optical lattice, interacting 
through a contact interaction: 
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d 3 r*t(i 



H + V conf +^(r)*(r) 



*(r), (1) 



where g = Airh as/m is the three dimensional interac- 
tion strength between bosons of mass m and scatter- 
ing length as- The single particle part of the Hamil- 
tonian defines the motional degrees of freedom through: 
Ho = — j^-y 2 + ^iat( r )' an d the confinement poten- 
tial: V con f = m[uj 2 (x 2 + y 2 ) + u> 2 z 2 ]/2. V CO nf defines 
the dimension, d, of the system. For d = 1 or 2 
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FIG. 1: The ratio of interaction matrix elements versus lat- 
tice depth in the lowest (dotted line) and the third (solid line) 
one dimensional bands. The inset shows the square of the 
Wannier functions in the lowest (dotted line) and the third 
(solid line) bands plotted as a function of distance against a 
confining, sinusoidal lattice. 



we have uj^-d <C u>d- The optical lattice potential is 
V£ t (z) = Vl [1 — cos (2irz/a)] /2, for d = 1. Here the 
lattice constant is a = A/ 2, where A is the wavelength of 
the laser defining the lattice. We also consider a square 
lattice for d = 2: V&(x,y) = V^{x) + V^{y). With 
these single particle potentials, the noninteracting prob- 
lem separates. In the direction of strong confinement 
we, as a first approximation, assume the harmonic oscil- 
lator ground state, separated from higher energy levels 
by Hud thereby establishing a d-dimensional problem in 
the remaining coordinates. The dimensionless interac- 
tion strength becomes ga — (g/Efi)(muJd/h)^~ (n/a) d . 
The over- bar indicates dimensionless units: f — irr/a 
and Ho = H /Er, where Er = h 2 /2m\ 2 . 

Along the directions of weak confinement the non- 
interacting problem defines a Bloch equation (exclud- 
ing the confinement potential). The exact solutions, 
&k.a, can be written in terms of Mathieu functions 
with wavevector k in band a 12]. From the Bloch 
functions we define the Wannier functions: Wi tOL — 

N» 5Z fe exp (—zfe • 5j)$fc. a (f), for AT S sites at locations 
Si. The Wannier functions localize in the "atomic" limit 

for large lattice depths, (Vl/Er)z — > oo, where the 
bands reduce to the harmonic oscillator energy levels. 
In a band, energetically, near the lattice maximum the 
density of two Wannier functions in neighboring sites can 
have strong overlap. The inset of Fig. ^ depicts two sit- 
uations, showing the square of the Wannier functions for 
d = 1 in the bands a — (dotted line) and a = 2 (solid 
line), plotted as a function of distance against a host 
lattice with height Vl — 2QEr. The large overlap be- 
tween nearest neighbor basis states suggests that atoms 
confined to higher bands offer the unique possibility of 
generating extended range interactions from an underly- 
ing, short-range interaction. In what follows we apply 



this technique to construct an extended, bosonic lattice 
model in the Wannier basis. A recent proposal [l3j sug- 
gests that extended range, Bose-Hubbard models may, 
alternatively, be generated with condensates of dipolar 
bosons in optical lattices. 

We now expand the operators in Eq. JIJ via: 'F = 
J2i a u>i,abi,a, where bi_ a annihilates a boson at site i in 
band a. The Hamiltonian becomes: H w = J2 a H a + 
J2 a =ia' H a a . We first focus on the largest intra-band 
terms: 

II" £ (4 A« + h - c ) £ n i>a 

<i,j> i 

+U a ' a ]T n i<a (n itCt - 1) + 4 V% a ni, anj , a , (2) 

i <i,j> 

where the chemical potential is defined to be /i a = 
Mo — (wi,a\Ho\wi t a) , with /l«o a constant dependent on 
confinement. The number operator is given by: rii jCt = 
b\ a b ia . The hopping between nearest neighbors, de- 
noted < i,j >, is only intra-band and non-diagonal for 
cubic lattices: tfj = — (wi <a \Ho\wj ja ). It is in princi- 
ple renormalized by the interaction to include conditional 
hopping: tfj -> tfj - 2Mf^ i f a (n l . a + n. ha - 1), where 
M ?i 1 ,QsTi' a * = (9dE R /2)(w iuai ;w l2 , a2 \w l3 , a3 - 1 w i ^ ai ) 
are the interaction matrix elements. In our study, we 
concentrate on the low density regime, p < 1, where p is 
the average number of particles per site. In this regime 
we find the conditional hopping to not change our results 
significantly. Along these lines we have, in Eq. (J2J), omit- 
ted double occupancy terms of the form &t a fet a b i a b i a 
which, as we have also checked, do not contribute sig- 
nificantly at low densities. The remaining two terms in 
Eq. J2J define the largest contributions to the interaction 
through the on-site, U a ' a = M^f c i ' i oc ' oc , and nearest 

neighbor, V°^- a = M° c '° L i f c ' <x , coefficients. Fig. ^plots 
the ratio V" i 'f' 1 /U a ' a as a function of the lattice depth 
for the lowest (dotted line) and the third (solid line) band 
in one dimension. In d = 2 the result remains the same 
as long as we compare the (0, 0) and (2, 2) bands. From 
Fig. ^ we clearly see that the ratio can be sizable. We 
must therefore incorporate extended Hubbard terms into 
any lattice model characterizing particles in higher bands 
not in the atomic limit. 

We now discuss a four stage gedanken experiment de- 
signed to place bosons in a higher band of the optical 
lattice. The prescription we provide here is not unique 
but serves to minimize band mixing. We first consider a 
partially filled lowest band in the atomic limit and with 
weak interaction strength, g C 1, As we know, from 
the mean field phase diagram 01 the Bose-Hubbard 
model (V — in Eq. J5J), bosons, in this limit, form a 
superfluid at alH/C/. 

The second step consists of adiabatically loading [f| 
the atoms into a higher band, e.g. a p — (2, 2) in d = 2, 



by oscillating the lattice depth at a frequency matching 
the inter-band energy difference. We assume that a large 
majority of the atoms can be transferred from the lowest 
band to a single, higher band. Once loaded into a higher 
band we note that, in our model, there is no inter-band 
coupling for a translationally invariant, non- interacting 
system in the steady state. 

In the third stage the lattice depth lowers, away from 
the atomic limit, to a point where the Wannier func- 
tions have some extension into the barriers between sites, 
Vl ~ 19-Er in Fig. ^ This process may be considered 
adiabatic if the time scale associated with lowering the 
lattice depth of a noninteracting system is much longer 
than h/\fj, a — p ap \, where a indicates the nearest band. 
We find \p a — /i( 2 ' 2 )| to cross zero linearly as a function 
of V L near V L w 15.8E R and V L « 17.5E R for a = (3, 1) 
and (4,0), respectively 

In the last stage we increase g~d- Recent studies 
0,0,0 have, quite differently, reached the strongly cor- 
related regime, in the lowest band, by tuning the ratio 
t/U with Vl- We, however, require the lattice depth to 
remain in a narrow regime. We assume that the interac- 
tion strength itself can be tuned through, for example, a 
Feshbach resonance. In what follows we study Eq. J5J, 
in the range g~d ~ 100 — 0. We then analyze inter-band 
effects induced by large interaction strengths. 

We consider several possible ground states of Eq. i ffl l 
and focus on the two dimensional square lattice. The 
ground states of this model contain four types of or- 
der, in the absence of disorder and at zero temperature: 
Superfluid order < bi >, Mott order < n, >, checker- 
board density wave order (— ) [< m > — p] , and 
supersolid order, where superfluid and density wave or- 
der coexist. Nonzero superfluid order arises from a spon- 
taneously broken gauge symmetry. We note that the 
host lattice corrugates the superfluid density at wavevec- 
tors corresponding to the reciprocal lattice vector. When 
phase fluctuations become strong Mott order persists at 
integer densities. 

The extended interaction term frustrates the Mott and 
superfluid phases leading to spontaneous translational 
symmetry breaking, e.g. the p = 1/2 density wave phase, 
ordered at half the reciprocal lattice vector. In contrast, 
a deep host lattice induces the corrugation in the Mott 
and superfluid phases. The fourth phase, the supersolid, 
arises from dual spontaneous symmetry breaking (both 
translational and gauge symmetry) inherent in coexisting 
density wave and superfluid order. 

We now discuss our solution of Eq. [21 in the band 
otp = (2,2). We solve Hq exactly to obtain the ma- 
trix elements. We use a Gutzwiller variational ansatz 
E| equivalent to a mean field decoupling of H a »: 



where the variational 



E^ i( „=o /jv<,„ \ N i,a) 
parameters, /jv 4 a , may vary over distinct sublattices 
and weight Fock states with N% a particles. We mini- 
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FIG. 2: Zero temperature, mean field phase diagram of the 
extended Bose-Hubbard model, Eq. J5J , as determined by the 
band structure of a two-dimensional, square, lattice in the 
(2, 2) band with interaction strength g2 = 50. 



mize Eq. (0) with respect to a keeping enough Ni. a 
to ensure convergence of the total energy. Note that 
Uj,p, Vi,j, and U depend only on V L /E R , g d , and po/E R . 
Fig.|21shows the two dimensional phase diagram for p < 1 
in the principal band, (2,2), with §2 = 50. The y-axis 
plots the chemical potential and the x-axis the lattice 
depth, both in units of E R . The hopping and the ex- 
tended Hubbard coefficients decrease with increasing lat- 
tice depth. Accordingly, we find Mott order at large lat- 
tice depths. Supersolid and density wave order appear 
f° r 92 ^ 40. The supersolid phase appears upon dop- 
ing of the density wave phase at p = 1/2 and not the 
Mott phase, consistent with the results of Ref. We 
add that in one dimension the phase diagram is nearly 
identical for the same set of parameters. However, it is 
by now well established that fluctuations destroy su- 
persolid order in one dimension. In both one and two 
dimensions (though more so in two dimensions) nearby 
bands energetically approach the principal band at low 
lattice depths. 

We now study interaction induced, inter-band effects. 
Our single band approximation, Eq. (0) , comes into ques- 
tion as we lower the lattice depth. We study, for d = 2, 
mixing with the two nearest bands (3,1) and (1,3). In 
principle, mixing with nearby bands can alter the phase 
diagram. However, if only a small fraction of the atoms 
occupy neighboring bands we may then safely assume 
that the phase diagram remains qualitatively the same. 
We ask whether or not the ground states of Fig. in the 
principal band at p = (2, 2 ) fa vor scattering processes cou- 
pling neighboring bands |l9| . The dominant inter- band 
terms are: 
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FIG. 3: The maximum probability that two particles from 
the (2, 2) band occupy the same site and scatter to the (3, 1) 
and (1,3) bands versus lattice depth with p = 2.2Er and 
g 2 = 50. 



where the matrix element in the last term ensures con- 
servation of band index (arising from conservation of lat- 
tice momentum). We have A = (—2,2) for a, a' 6 
{(2, 2), (1, 3), (3, 1)} in two dimensions. 

The last term in Eq. (PJ) takes two particles from 
the same site in the principal band and "scatters" 
them to neighboring bands and vice versa, when ap- 
plied to a state initially in a p . It acts as the domi- 
nant inter-band scattering mechanism |20j . We calcu- 
late the probability of such an event through first or- 
der perturbation theory. Consider two states ijj ap and: 

= ■^- 1 (^ ap +A/2 & U p -A/2^, ap ^, ap )^« p , where M 
is a normalization constant. In the absence of dissi- 
pation, the probability that two particles at any one 
site occupy neighboring bands oscillates in time, T: 
P(T) = 2\{^ e \H w \4, ap )\ 2 [I - cos (ET/h)}/E 2 , where 
E = {np e \H w \tj) e ) - (ip ap \H w \ip ap ). We argue that if the 
probability remains small, then band mixing will be sup- 
pressed. Note that, with a large interaction strength, E 
is not equal to the single-particle, self energy difference 
between bands. 

The probability of finding two particles at the same 
site is small in all regions of Fig- El Ignoring fluctuations, 
the p — 1 Mott and p = 1/2 density wave phases have 
no double occupancy. Therefore, the superfluid and the 
supersolid phases remain as the only candidate phases 
involving on-site, inter-band scattering processes. Fig. [31 
plots the maximum scattering probability at T — h/2E 
as a function of lattice depth for a chemical potential 
M (2,2) = 2.2E R and g 2 = 50. The probability of finding 
two particles in nearby bands is less than 2 x 10~ 2 . This 
suggests that at intermediate lattice depths and low den- 
sities the single-band ground states resist on-site scatter- 
ing processes into neighboring bands. The energy differ- 
ence, E, remained non-zero for all Vl primarily because 
of the strong inter-band, nearest neighbor interaction in 
the neighboring bands. 

We have shown that promoting bosons to higher bands 
of optical lattices can lead to states beyond the superfluid 
and Mott states present in zero-range, Bose-Hubbard 



models of the lowest band. We argue that the ground 
states of an extended Bose-Hubbard model capture the 
essential physics of bosonic atoms placed in a single, 
isotropic band with a minimum lying, energetically, near 
the top of the optical lattice. The resulting supersolid 
and density wave states add to the set of observable 
phases. 

Additional phases may arise outside the set of approx- 
imations leading to Fig.|3 Gutzwiller, mean-field theory 
should be an excellent approximation for d > 1 and af- 
firms results obtained from quantum Monte Carlo stud- 
ies for d = 2 21] . However, our results overestimate 
the strength of the supersolid phase because we have ex- 
cluded a competing phase-separated state j^. Further- 
more, strong inter-band mixing can populate anisotropic 
neighboring bands (e.g. (3,1) and (1,3)) leading to stripe- 
like superfluid states which coexist with the superfluid 
and supersolid states in the primary band. 

The states confined to band a p are, technically, 
metastable. We require r 3> h/\t i ^ p \ 1 where t is the life- 
time of the state, r may be affected by dissipative effects 
including collective mode inter-band scattering. Our re- 
sults for the realizable extended Bose-Hubbard Hamil- 
tonian, in addition to providing a rich quantum phase 
diagram, yields an interesting connection to topological 
quantum computation Q in cold atom optical lattices. 
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